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We analyze the statistics of the electromagnetic radiation emitted from electrons pushed through a 
quantum point contact. We consider a setup implemented in a two-dimenional electron gas (2DEG) 
where the radiation manifests itself in terms of 2D plasmons emitted from electrons scattered at the 
point contact. The bosonic statistics of the plasmons competes with the fermionic statistics of the 
electrons; as a result, the quantum point contact emits non-classical radiation with a statistics which 
can be tuned from bunching to anti-bunching by changing the driving voltage. Our perturbative 
calculation of the irreducible two-plasmon probability correlator provides us with information on 
the statistical nature of the emitted plasmons and on the underlying electronic current flow. 



I. INTRODUCTION 

The interest in photon radiation from a quantum point 
contact (QPC) is two-fold: on the one hand, the quantum 
point contact acts as a source of non-classical lights. 
While a classical current produces photons with Pois- 
son statistics^, the current pushed through a QPC can 
be tuned to radiate photons with either super- or sub- 
poissonian statistics^, testifying for the bosonic na- 
ture of the photons (bunching) or the fermionic statis- 
tics (anti-bunching) of the underlying electrons radiat- 
ing these photons, e.g., as 'Bremsstrahhmg' radiation 
in back-reflection processes. On the other hand, the 
radiation statistics contains information on the statis- 
tics of current fluctuations across the QPC — as such, 
the photodetector serves as a tool to probe the nearby 
current flow in the mesoscopic wire^^. Both themes 
have attracted considerable interest, be it in the context 
of the counting statistics of electrons in phase coherent 
mesoscopic conductors'^, where higher-order correla- 
tors carry the signatures of interactions, or be it related 
to the search for new sources of non-classical radiation 1 ' 2 . 

The classic photodetection theory for optical photons 
goes back to Glauber— and is based on a threshold detec- 
tor, where each (sufficiently energized) photon induces an 
electronic cascade generating a counting signal. The the- 
ory has been widely applied in the optical regime where 
threshold energies reside in the eV regime. The anal- 
ysis of the photon counting statistics in Refs. has 
been based on this photodetection theory and has con- 
centrated on the regime, where the energy of emitted 
photons is larger then eV/2, with V the voltage bias ap- 
plied to the quantum point contact. For this situation, 
the emitted photons are anti-bunched, since they are pro- 
duced by different electrons scattered by a QPC. 

The application of Glauber's theory to mesoscopic se- 
tups is not straightforward, however, as typical photons 
generated in a mesoscopic structure have energies in the 
GHz regime and it may be difficult to construct a sta- 
ble photodetector with a correspondingly small threshold 
energy. Furthermore, the emission rate of free photons 
scales with the small factor a(v F /c) 2 ~ 10~ 8 , where a is 
the fine structure constant, v F is the Fermi velocity, and 
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FIG. 1: Quantum point contact (QPC) fabricated through 
gates constricting the electron flow in a two-dimensional elec- 
tron gas (2DEG). We consider a quasi- ID setup where the 
width w of the 2DEG allows propagation of one transverse 
mode; d is the distance to the backgate. The voltage (V) 
driven quantum point contact emits 2D plasmons with fre- 
quencies fii and fi2 due to electron backscattering (t and r 
denote transmission and reflection amplitudes). These plas- 
mons are picked up by waveguides transmitting the signal to 
the measurement setup. The measurement of the one- [Pi (fi)] 
and two-plasmon [Pz(fli, 0.%)] emission probabilities allows to 
characterize the radiation and provides information on the 
fourth-order current correlator. 



c the velocity of light, rendering the observation of free 
photons practically impossible. 

On the other hand, within a mesoscopic context, the 
low-energy exchange between a quantum point contact 
and the electromagnetic environment can be studied with 
a double-dot detector^i^. Furthermore, a setup has been 
analyzed recently, where low-energy radiation manifests 
itself in terms of plasmons propagating within the two- 
dimensional electron gas (2DEG) forming the quantum 
point contac t 11 ! 12 , cf. Fig. Q] These 2D plasmons then 
are further transmitted through waveguides and then an- 
alyzed in a power detector. The emission rate for plas- 
mons is enhanced over that of free photons by a factor 
(c/upi) 3 ~ 10 6 , where v p \ denotes the plasmon veloc- 
ity; this gain in signal has to be preserved by proper 
impedance matching of the waveguides and the quan- 
tum point contact. In this paper, we study this kind 
of setup and analyze the two-plasmon correlations in an 
arbitrary frequency range; these correlations then tell 
about the statistics of the emitted radiation (bunching 



2 



or anti-bunching) and carry information on the fourth- 
order current correlator in the mesoscopic wire. 

Below, we present a perturbative calculation (in the 
QPC-electromagnetic field coupling) of the one- [Pi(fi)] 
and two-plasmon [P2(fii,f22)] probabilities for emission 
at given frequencies f2, O12 and during an arbitrary but 
fixed time to- Within our perturbative approach, we cal- 
culate the irreducible probability correlator p^i, ^2) = 
P 2 (fti,fi 2 ) - Pi(^i)Pi(^ 2 )- The sign of this quantity 
then tells us about the character of the emitted radiation, 
bunching or anti-bunching. E.g., a positive sign indi- 
cates that plasmons prefer to be emitted simultaneously 
rather then independently and thus are bunched, while 
for a negative sign the opposite situation of anti-bunching 
prevails; for P2(f2i,f22) = the two plasmons are emit- 
ted independently. As we will show below, changing 
the bias on the quantum point contact will allow to 
change the statistics from anti-bunched fermionic-type 
to bunched bosonic-type. Furthermore, for different fre- 
quencies ill 7^ Q2, the quantity ^2(^1,^2) is directly 
proportional to the irreducible fourth-order current cor- 
relator and hence its measurement provides valuable in- 
sight into the fluctuation statistics of the current across 
the quantum point contact. 

In the following, we first define the model and find 
the expressions for the probability densities pi(f2) and 
732(^1,^2) expressed through current correlators within 
a perturbative expansion, see Sec. [TTJ In a second step, 
these probability densities are rewritten through the sec- 
ond and fourth order irreducible noise correlators 
and S( 4 \ In section ITTT1 we combine results for the prob- 
ability densities and the noise correlators to find the sin- 
gle photon emission probability Pi(f2) = px(Cl)dfl, the 
correlated two-photon emission probability P2(f2i, f2 2 ) = 
p 2 (yii, f^df^df^, as well as the irreducible quantity 
P^i,^); we identify the regimes of anti-bunching at 
intermediate voltages and the competition between anti- 
bunching and bunching at high voltage, when one elec- 
tron has sufficient energy to emit two (bunched) plas- 
mons. In section lTVl we discuss the interrelation between 
the fourth-order current correlator and the irreducible 
correlator ^2(^1,^2) and we conclude in Sec. IVl 



II. PHOTON EMISSION PROBABILITIES 

In this section, we derive the formal expressions for the 
probability densities pi(f2) and p 2 (fii, ^2) to emit one 
or two plasmons at frequencies f2, Qi i2 during a given 
time and express the results through the second- and 
fourth-order current correlators. We restrict ourselves to 
a perturbative analysis in the Hamiltonian Pint, describ- 
ing the interaction between the electronic current density 
j(r) through the QPC and the electromagnetic field de- 
scribed through the vector potential A(r), 

Hint = -\j d 3 r j(r) • A(r). (1) 



In a realistic setup, the electromagnetic field is gen- 
erated by the two-dimensional plasmon modes of the 
two-dimensional electron gas (2DEG) 13 . These plasmon 
modes propagate along the one-dimensional wire (the x 
direction) and we consider the lowest transverse mode 
only (we assume translation invariancc along the y direc- 
tion). The presence of the backgate changes the disper- 
sion from the usual (in 2D) square-root form to a linear 
one. The vector potential A.(x) then has only a compo- 
nent along x, 

A x {x) =^k 1 (^^-) 1/ \c k e^ - gt «-*■), ( 2 ) 

\LUk -Equant J 

where cl (ck) are bosonic creation (annihilation) opera- 
tors for the (longitudinal) plasmon modes with wave vec- 
tor fc; ojk — v p \k is the plasmon frequency with v p \ the 
plasmon wave velocity and 7 ~ \Jdjw is a geometrical 
factor with w the width of the 2DEG and d the distance 
between the 2DEG and a Schottky barrier electrostatic 
gatei 3 -. In the following, we set the quantization length 
Equant equal to unity. 

With the current operator j x (r) — I(x)8 5 (v — le^) in 
the wire, the Hamiltonian P; nt in the interaction repre- 
sentation takes the form 

Pint(i) = -*7$> ( 3 ) 

x(cj k {t)e-^ -tlility^), 

where Ik(t) is the spatial average of the current I(x,t) 
over the coupling region of the plasmon described by the 
kernel f(x) with extension L, 

I k {t) = J dxf(x,t)f(x)e ikx . (4) 

In a typical situation, the frequency of the excited plas- 
mons lies in the GHz range, with the velocity v p \ roughly 
100 times slower than the speed of light c. The corre- 
sponding wave length A p i ~ 100 /im then is much larger 
than the size ~ 0.1 p oi the quantum point contact. As- 
suming a coupling region 0.1 pm < L < A p i, we can 
ignore the k dependence in ifc(t). 

We assume that initially, at time t = —to < 0, no 
plasmons are excited and the state of the total system 
(QPC and bosonic modes) is described by the factorized 
density matrix p{—t ) — p p i(— to)®P qPC (— to)- At a later 
time t = 0, the probability density p n (f2i, . . . , fl n ) to find 
n plasmons with frequencies . . . , f2„ can be defined in 
terms of the time-ordered evolution operator S(0, —to) — 
Texp[-(i/H)J° to H int (t')dt'}, 

p„=Tr{p„(Q 1 ,...,n n ) ) §(0,-to)/3(0)5t(0,-to)} J (5) 

where pit) is the time dependent density matrix 
of the electronic system (including scattering at the 
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QPC, interactions between electrons, etc.), V n = 
\k\, . . . , k n ){ki, . . . , k n \ is the projection operator on the 
state with n bosons of frequencies tti = v p \ki, and the 
trace is taken with respect to the total system (note that 
V n {Q) = V n {k)Y\S((l — v p \k); the projection does not 
select specific directions of the emitted plasmons). The 
Taylor expansion of the evolution operator S(0, t) up to 
the lowest non-trivial order provides us with the (one 
plasmon) probability density, 



n 



Pl (fi) = 7 2 ^ J dsdre iQ ^ (I(s)i(r)), (6) 

P -to 

where (...) denotes the average over the electronic sys- 
tem. The next term in the series generates the two- 
plasmon probability density 
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J dsi ds 2 drx dr 2 (7) 



x(T_{f( J » 1 )/( fl2 )}r + {7(r 1 )/(75)}) 
x (e 



,tni(n — Sl) p in 2 (-T2-S2) I gt£Jl(n— S2) p <n 2 (T2 — Si) 



-f2i <-> f2 2 ), 



where T + and T_ are the time ordering operators in the 
forward and backward directions, respectively. In the 
following, we assume to 3> f^ 1 ,^^ and regularize the 
time integrals in Eqs. ([6]) and 10 at the lower limit by 
introducing a small damping factor exp(— 77|r|) with rj <C 
fL Physically, t]~ 1 corresponds to the decay time of the 
plasmon excitations propagating in the wave guide. 

The above time integrals over current correlators can 
be expressed through the spectral power of current fluc- 
tuations. Assuming a stationary situation, the current 
correlators in Eqs. ([6]) and {7} depend only on relative 
times; the second-order noise correlator then can be 
defined through 



dr e 



((i(r)I(0))), 



(8) 



and the fourth-order correlator reads 



(u>i , lo 2 , ) = / dri dr 2 drz e 



- 1U) 1 Ti — IUJ2 T2 — T3 



x «J(ti + t 2 + t 3 )I(t 2 + 7^)1(73)7(0))). (9) 

The noise correlators and involve irreducible 
current correlators, while the expressions for Pi(Q) and 
p 2 (fil, f2 2 ) are expressed in terms of reducible quantities. 
Expressing the reducible correlator through irreducible 
ones, 



(7( S )7(r)) = ((7( S )7(r))) + (7( S ))(7(r)), 



(10) 



and using the Fourier transform ((7(s)7(r))) = J(dui/2w) 
exp[iuj(s — t)]S^- 2 \u>), we can perform the time integrals 
in the expression Eq. Q for pi(Q). In the stationary 



regime, the average current (I(t)) = I through the QPC 
is independent of time and we obtain the intermediate 
form 



I 2 



dm S^{uj) 

2tt (w-17) 2 + 77 2 " r f2 2 +?7 2 



(11) 



In the limit r\ -C f2, we can approximate the Lorentzian 
by a (5-function, r/ / [(qj — ft) 2 + r] 2 ] ss itS(uj — £1), and carry 
out the integral over w to arrive at the final result, 
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(12) 



Furthermore, the contribution from the irreducible part 
of the current correlator to pi{fl) is a factor Sl/77 ^> 1 
larger than the contribution from the reducible part cx I 2 , 
as follows from the estimate (Q) ~ I 2 /Q we will drop 
the reducible part in our further analysis below. 

Expressing p 2 (Vti, VL 2 ) through irreducible correlators 
is more involved: below, we keep only those terms of 
the reducible fourth-order current correlator which give 
non-vanishing contributions to p 2 in the limit rj <C fi, 

(7( S2 )7( Sl )7(n)7(T 2 )) = ((7( S2 )7( Sl )7(ri)7(r 2 ))) 
+((/> a )/(r 2 )))((/( ai )/(r 1 ))) + ..., (13) 

where we assume t\ > r 2 , si > s 2 and . . . denotes terms 
with different time orderings as well as third-order cumu- 
lants providing irrelevant contributions. The probability 
distribution p 2 can then be expressed through the noise 
correlators and S^; the contribution from the re- 
ducible part of the fourth-order current correlator reads 



p i 2 2 \n 1 ,n 2 ) = Y 



X 1 + 



fii n 2 s(- 



>(ni)^ 2) (n 2 ) 



H H 

V 



Arj 2 



(fii - n 2 ) 2 + if 



(14) 



This term describes the excitation of two plasmons due 
to the uncorrelated scattering of independent electrons. 
For two bosons with equal frequency |fii — fi 2 | <C rj, 
this contribution is enhanced by a factor 2 as compared 
with the probability to emit two bosons with different 
frequencies (see Eq. ([T2"j)). 



Pi(fii)pi(n 2 ), |Qi-n 2 |>??, 
,2/0 \ ,0 o 1 ^ _ (1 5 ) 



2pf("i 



l^i - n 2 | < 



This enhancement is a quantum mechanical time- 
interference effect: for |f2i — f2 2 | -C rj we cannot know 
which boson was emitted first during the measurement 
time rj -1 and the amplitudes of both alternatives have to 
be added, resulting in a constructive interference between 
them. 

Next, we concentrate on the contribution p 2 arising 
from the irreducible part of the fourth-order current cor- 
relator in Eq. (|13[) . After integration over times, we arrive 
at the expression 
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. (wi — Oi — W7) (c^3 — f2i + iry) (a>i — fix — iff) ( w 3 — ^ 2 + «?7) 



fix fi ; 



(16) 



Combining Eqs. (|14[) and (JTHJ) and approximating the Lorentzian in Eq. (|16p by a (5-function, one finally arrives at 
the expression for the probability density to emit two plasmons, 



with 
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FIG. 2: Setup to probe the emission of photons from a QPC 
into free space. The photons are picked up by the LC cir- 
cuits 'up' and 'down', which are inductively (i u (d),qpc) cou- 
pled to the QPC and capacitively (C u (d),dot) coupled to the 
quantum dots. The latter are tuned to minimal or maximal 
transmission and thus operate as quadratic detectors. Pho- 
tons absorbed in the LC-circuits modulate the charge on the 
capacitors C u (d) and, through coupling to the quantum dots, 
the probe currents 7 U and J d . 



The resulting probability densities Eqs. (| 1 2|) and (TT7|) 
involve second-order current correlators at positive fre- 
quencies and integrals over fourth-order current correla- 
tors with integration kernels concentrated near positive 
frequencies as well. This feature is easily understood: ini- 
tially there are no bosonic excitations, hence within our 
lowest-order calculation the only processes contributing 
to Pi(fl) and p 2 (ili,i7 2 ) are due to plasmon emission — 
reabsorbtion processes involving negative frequency cor- 
relators show up only within a higher-order analysis. 

Finally, we note that our results hold true, up to a 
numerical factor of order unity, for the emission of (long 



wave length) photons as well; a corresponding setup with 
LC circuit pick-ups and quadratic detectors is sketched 
in figure O For the photonic emission, one needs to re- 
place the plasmon wave velocity v p \ by the speed of light 
c and set the geometry factor 7 = 1 in Eqs. (fT2"]) and 
(jTTJ) . We will see below, cf. Eq. [2H that the basic dimcn- 
sionless parameter governing the numerical value of our 
results is a cx (e 2 /hv p i)(vF/v p i) 2 . As a consequence, the 
one- and two-photon emission rates are suppressed by the 
small factors (u p i/c) 3 (with v p \/c ~ 10 -2 typically) and 
(vpi/c) 6 , respectively, when real three-dimensional pho- 
tons are emitted into space. In the following, we use the 
terms plasmons and photons synonymously. 



III. PHOTON COUNTING STATISTICS FOR A 
QUANTUM CURRENT 

It is a well known result due to Glauber 3 that a classi- 
cal current I(t) produces a coherent state of the electro- 
magnetic field, with the width of the Poisson statistics 
of photo counts for each mode given by the correspond- 
ing Fourier coefficient I(lo) of the current. This implies 
that photons are emitted independently, i.e., the joint 
probability to emit two photons with frequencies f2i and 
i7 2 is in fact a product of single photon probabilities, 

p 2 (fii,n 2 ) = Pi(fii)Pi(n 2 ). 

For a quantum current the fluctuations appear due to 
the scattering of separate electrons. In the scattering pro- 
cess, the electrons experience an acceleration and emit 
Bremsstrahlung radiation. The fermionic statistics in- 
duces a correlated flow of the electrons incoming from the 
voltage biased reservoir and separate electrons are scat- 
tered one by one. The photons emitted in the scattering 
of separate electrons then inherit their fermionic correla- 
tions. On the other hand, one electron may produce sev- 
eral photons during the scattering event, in which case 
these photons are bunched. The deviation from poisso- 
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nian statistics thus is a competition between these two 
processes. 

Using the results of the previous section, we can find 
the probability Pi(fi) = pi(f2)df2 for the emission of 
plasmons in the frequency interval d£l around O and 
f^(Oi, fi 2 ) — P2(fii, ^2) driidil2 for the emission of two 
plasmons at given frequencies within the non-overlapping 
frequency windows dfli and e?f2 2 around fij, f2 2 ; for over- 
lapping frequency intervals Q 1 « J7 2 — the determina- 
tion of P 2 (f^) requires proper integration of the density 

p 2 (Oi,fi 2 )- 

While the results on the one- and two-plasmon pro- 
cesses provide us with only limited information on the 
photo counting statistics, they are nevertheless sufficient 
to tell us about important quantum signatures in the ra- 
diation. In particular, the statistical properties of the 
emitted plasmons is conveniently described by the irre- 
ducible probability correlator, 

P 2 (Q u n 2 ) = p 2 (n u n 2 ) - p 1 (o 1 )p 1 (o 2 ). (19) 

A positive sign of P 2 (Qi, f2 2 ) indicates that plasmons are 
bunching, i.e., they are preferentially emitted simulta- 
neously. In the opposite case of negative correlations, 
P 2 (£!i, Q 2 ) < 0, the plasmons are anti-bunched, implying 
that the probability to find the second plasmon emitted 
right after the first is suppressed. 

A. Single plasmon emission 

In the following, we assume that the emitted radiation 
is collected from a finite region of length L < X p \ and we 
average the current operators in Eq. (fjQ) using a specific 
kernel of the form f(x) — exp(— \x\/L). First, we concen- 
trate on the single plasmon probability Pi(il). Inserting 
the (T = 0) second-order noise correlator- 
s'^) = —TRe l(n/v ^ X2 ^ x ^{uj -VL)Q{Lo Q -VL) 

(20) 

into Eq. (TT2"]) (Q(x) is the Heavyside function), one ar- 
rives at the result for the single plasmon emission prob- 
ability in the form, 

(21) 

with u>o — eV/h the voltage frequency, A = v-p/Q ~ 
100 jj,m the characteristic spatial scale of the current fluc- 



tuations, and a is the dimensionless parameter, 




Note that with d and w of order 1 /xrr>i£, hence the ge- 
ometry factor 7 is of order unity. On the other hand, the 
parameter a is reduced due to the residual impedance 
mismatch between the quantum point contact and the 
waveguides, cf. Fig. [TJ 

The probability P\(Vl) describes the Bremsstrahlung 
plasmon emission due to electron scattering at the bar- 
rier. Although one may expect the probability Pi(f2) to 
increase for a more effective scatterer with T — > 0, the re- 
sult is in fact proportional to T(\ — T) and vanishes in the 
tunneling limit. This peculiarity of the Bremsstrahlung 
appears due to the Fermi statistics: the electron com- 
ing in from the biased reservoir has to relax to a state 
with a lower energy after the plasmon emission. Since 
all electron scattering states of the biased reservoir with 
lower energy are filled, the only possibility for the elec- 
tron to decay is into an empty scattering state of the un- 
biased reservoir. This process requires tunneling of the 
electron through the barrier and hence the probability 
Pi(f2) turns out proportional to T. The sharp suppres- 
sion of Pi(f2) when the plasmon energy Ml is larger than 
the applied bias eV can be explained in the same way: 
an electron cannot find an empty state to emit such a 
'high-energy' plasmon. 



B. Correlated plasmon emission 

Next, we analyze the two-plasmon probability 
P 2 (f2i,f2 2 ), see Eq. (fT7|) . which involves irreducible 
fourth-order noise correlators. Below, we consider the 
case of an extended interaction region, where the signal 
is collected over a region of size L larger than the char- 
acteristic spacial scale A = uf/^i,2 ~ 100 /im of current 
fluctuations, where we consider typical Fermi velocities 
vp ~ 10 7 cm/s in GaAs heterostructures and frequen- 
cies in the GHz regime. The calculation of the expres- 
sion (Qi, f2 2 ) in Eq. (|18|) is carried out in the Ap- 
pendix, cf. Eq. (|A19[) for the final result; the probability 
P 2 (f2i,f2 2 ) for two-plasmon emission in non-overlapping 
frequency intervals then reads 



p 2 (n 1 ,n 2 ) = —{RT(T~R) 2 e(La -n s )(Lo ~ni:)-2(RT) 2 Q(Lj ~n > )(Loo-n > ))^^ 



(23) 



4q 2 
(tti])'- 



(PT) 2 e(w - Hi)e(a;o - J2 2 )(w - «i)(^o - ^2) 



1 



dfli c?f2 2 
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where Sis = Sli + VI2 and f2> = maxjili, ^2}- The first term oc I/77 in the above expression originates from the 
irreducible fourth-order noise correlator while the second term oc I/77 2 is a contribution from the reducible part 

p!p of the fourth-order current correlator, see Eq. (fl~7|) . 

For non-overlapping frequency bands with |fii — O2 1 ^> 77, the second term in the bracket [. . .] in the reducible 
part is small and what remains is equal to the product Pi(f2i)Pi(£l2) of single plasmon probabilities. Hence, the 
irreducible probability correlator P 2 , Eq. (fl"9|) involves only the first irreducible term in Eq. ([23]) . 



irrj L 



dQidQ2 

^1^2 ' 



(24) 



In the opposite situation where the frequency bands 
overlap, the density ^2(^1,^2) has to be properly inte- 
grated over frequencies; the second term (of Lorentzian 
shape) in p 2 2 ^ , cf . Eq. (TT4")) , then provides a contribution 
tt/t]. For a narrow frequency band 5fl/fl <C 1 far from the 
voltage frequency \loo — Q\ ~ Q, this reducible contribu- 
tion dominates over the irreducible one originating from 
Q*- 4 -* and the irreducible probability correlator Eq. (fl~9]) 
is dominated by the reducible part of the fourth-order 
current correlator, 

P 2 (Q) « ^-(Ti?) 2 e(c - Sl)(uo - n) ^_R 6 J±. (25) 

TTT] il II 

The situation changes if one detects plasmon in the 
overlapping frequency band reaching the voltage fre- 
quency, f2i,f22 £ [u>o — Af2, wo], Af2 < luq/2. In this 
situation one needs to take the remaining integrals over 
the frequencies f2i,Q 2 in Eq. ([23]) exactly. The result 
takes the form 

P 2 (0) = -- (PT) 2 + — (RTf^r, (26) 

where the first term originates from the irreducible part 
of the fourth-order current correlator, while the second 
term has its origin in the reducible contribution. We then 
find that the irreducible contribution dominates in this 
regime. 

In summary, the irreducible correlator P2 at differ- 
ent frequencies f^i ^ f2 2 involves a competition between 
bunching (first term in (|24p ) and anti-bunching (second 
term in (f2~4"]) ). When the frequency intervals strongly 
overlap we encounter two regimes: 1) detecting the plas- 
mon in a narrow frequency band far from the voltage fre- 
quency, \loo — fl\ ^ fl ^ Sfl, the reducible contribution 
of the current correlator always dominates thus resulting 
in bunched plasmon radiation; 2) measuring the plasmon 
in the frequency band [ujq — AO, loq], Afl < wq/2, near 
the voltage frequency, the irreducible contribution of the 
current correlator for a single-channel conductor is two 
times larger then the reducible one in absolute values, 
resulting in anti-bunched radiation. 

It is the irreducible correlator P 2 at different frequen- 
cies f2i ^ f?2 which is the most interesting quantity mea- 
suring the non-trivial correlator Q^ 4 \ Thus, in order to 
study the fourth-order electron noise effects on the plas- 



mon statistics, one needs to be able to distinguish be- 
tween the frequencies of the emitted plasmon during the 
measurement time n^ 1 . 

We thus concentrate on the detection of plasmon with 
distinguishable frequencies. The probability correlator 
P 2 (Oi, fi 2 ), see Eq. (|2"4"]) . involves two terms with oppo- 
site signs, a positive one proportional to O(wo — fii — 2 ) 
and a negative contribution oc 0(luq — max{f2i, f2 2 }). Ap- 
plying a small voltage bias loq < max{f2i, O2}, the scat- 
tered electrons do not have the possibility to emit two 
plasmons with frequencies fli and O2 due to the restric- 
tion of the Fermi statistics on the final electron state and 
hence P 2 (Oi,0 2 ) = 0. 

In the intermediate voltage regime max{f2i,f22} < 
Loo < Oi + Q2, a single electron can emit only one plas- 
mon, either of frequency fii or f2 2 . Thus the corre- 
sponding contribution oc Q(loq — max{J7i,02}) is due 
to the scattering of different (fermionically correlated) 
electrons. At zero temperature, these electrons arrive 
at the scatterer in perfect order with a time separa- 
tion Ty ~ h/eV. Hence, the second plasmon is emit- 
ted only after the time Ty, resulting in an anti-bunched 
radiation statistics and a negative probability correlator 
P 2 (Oi,0 2 ) < 0. When a single electron creates only 
one plasmon, the complexity of the emission process is 
reduced considerably, allowing for the determination of 
the full counting statistics of the emitted radiation^. 

Increasing the voltage beyond h(Qi + ^2)7 the above 
single plasmon production is augmented by processes 
where one electron emits two plasmons in a single scat- 
tering event; this process generates bunched radiation 
and hence the corresponding term oc 0(luq — f2i — f2 2 ) 
is positive, see Eq. (f2~4"|). The overall sign of the probabil- 
ity correlator P 2 (Oi,f2 2 ) then appears in a competition 
between the two processes creating plasmons individu- 
ally or in pairs. Changing the parameters of the quan- 
tum point contact, one can control the relative weight 
of the two contributions. For example, at T = 1/2 the 
two-plasmon process oc TP(T — P) 2 vanishes while the 
single-plasmon term oc 2 (TP) 2 is maximal, resulting in 
anti-bunched radiation in the whole two-plasmon voltage 
regime fii + f2 2 < < ^1 + ^2 + minjfii, f^} (we con- 
sider three-plasmon processes involving at least one plas- 
mon with frequency Qi and one with f2 2 ). Alternatively 
in the tunneling limit T « 1, the emission of bunched 
plasmons is the dominant process at u>q > fl\ + fl 2 - A 
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further increase of the bias voltage may lead to multi- 
plasmon production processes, where more than two plas- 
mons are created by a single scattered electron. However, 
such processes appear only in higher orders of pertur- 
bation theory, while we have restricted ourselves to the 
fourth order in the coupling constant. 

Finally, we analyze the case of a multi-channel con- 
ductor in the diffusive regime. The probability correlator 
P 2 (9i,9 2 ) for non-overlapping frequencies in the multi- 
channel case can be obtained by summing the terms in 
Eq. (|24[) over the channel index n with appropriate values 
for the scattering coefficients T n and R n . The distribu- 
tion of the transmission eigenvalues T in the diffusive 
conductor is given by the Dorokhov distribution func- 
tion^, 



m = — 1 

P{ ' 2G TVT^T 1 



(27) 



where G is the average conductance and Go = 2e 2 //i. 
Averaging Eq. (|2~4")) over T, the probability correlator 
P 2 (9i,9 2 ) for the radiation emitted from the diffusive 
conductor is given by the expression (we assume !![ / 
"2), 

p 2 (9 1 ,9 2 ) = — £- (ii e(c.o - (a*, - 9 S ) 

TTT] Gq V 105 



4 
35 



9(o; -0 > )(a;o-0 > ))^ 



dQi d9 2 

TV 



(28) 



Quite remarkably, the result is negative, indicating anti- 
bunched plasmons, even in the two-plasmon regime luq > 
ill + ^2- This provides, besides the measurements of the 
average current and noise, another consistency test for 
the Dorokhov distribution function. 



C. Measurement setup 

Next, we relate the probabilities Pi (9) and P 2 (f2i, ^2) 
to physically measurable quantities. In a realistic experi- 
ment, see Refs. lTllll2l . a two-terminal quantum point con- 
tact is realized in a two-dimensional electron gas inserted 
between two transmission lines, cf. Fig. [TJ The excited 
plasmon excitations in the 2DEG induce an ac-electric 
current signal in the two (left and right) transmission 
lines. The transmitted signals are independently ampli- 
fied by two cryogenic amplifiers and then passed through 
frequency filters (selecting proper frequencies f2i, 2 ) fol- 
lowed by quadratic detectors. Thus the signal taken at 
the end of each transmission line is proportional to the 
power emitted from the QPC, 



Wi = A i h9 i h(9i)d9 l + Wi, i = l,2, 



(29) 



where Ai is the amplification factor, h(9)d9 is the plas- 
mon production rate, i.e., the number of excited plas- 
mons per unit time within the frequency band <if2, and 
Wi is a parasitic power due to the intrinsic noise of the 



amplifiers and the vacuum fluctuations of the bosonic 
mode. Given the life time (or detection time) 1/77 of the 
plasmons, cf. Eqs. (TT2")) and (TTTj) , the plasmon production 
rate n(9)d9 relates to the plasmon occupation number 
N(9) via n(9)d9 = r]N(9). 

Since the noise signals u>i j2 are n °t correlated for 
the different transmission lines, {5w\5w 2 ) = 0, dwi = 
Wi — {wi), the irreducible cross-correlator between the 
two transmission lines, Q12 = (W1W2) — (Wi)(W 2 ), in- 
volves only the irreducible cross-correlator of the plasmon 
occupation numbers emitted into each line, 



Q12 = A X A 2 ffijffla r, 2 «JV(ni)JV(n 2 )». 



(30) 



Next, we express the plasmon number correlator 
((jV(fti)jV(fi 2 ))) through the probabilities Pi(fi) and 
P 2 (9i,9 2 ) calculated above. Up to fourth-order in per- 
turbation theory, we have 

(N(fi)) ^ Pl (9)d9 + Jdn' P2 {n,n') <m, (31) 

(N(9 1 )N{9 2 )) MjtoCfii.ffcJdfixdfia, (32) 

where the second term in {JV(S1)) is a higher order cor- 
rection. As a result, the irreducible cross correlator for 
the plasmon number then assumes the simple form, 

((N(9 1 )N(9 2 ))) =P 2 (fi 1 ,Q 2 )-P 1 (0 1 )P 1 (fi 2 ) ) (33) 

involving only the probabilities P 2 (f2i,f2 2 ) and Pi(O); 
other terms are of higher order in the coupling constant. 
The cross correlator of the emitted power between two 
transmission lines then is proportional to the probability 
correlator P 2 (f2i,fi 2 ), see Eq. (TT9")) , and the sign of the 
power cross correlator directly characterizes the statistics 
of the emitted radiation. 

Finally, we connect our results with the results of 
Beenakker and Schomerus^, which are based on the 
Glauber formula for photon counting. Within this ap- 
proach the relevant physical quantity to observe devia- 
tions from the poissonian statistics is the variance of the 
detected particles: Var(A) = (N 2 ) - (N) 2 . Within low- 
est order perturbation theory, Var(iV) can be expressed 
through the probabilities Pi and P 2 , 



Var(iV) = 4P 2 + P!-Pi 2 
= (N)+2P 2 -P 2 . 



(34) 
(35) 



The sign of the combination 2P2 — P 2 quantifies the devi- 
ation of the photon statistics from the poissonian result 
Vai(N) = (N). Let us first concentrate on the regime 
where photons are measured in a narrow frequency band 
[9 - 59,/ 2, 9 + 59/2], 59 < fl, far from the voltage fre- 
quency \u)q — 9\ ~ 9. According to the results of Ref. U 
the quantity 2P2 — P 2 is proportional to the measurement 
time t = I/77, i.e., terms oc l/?7 2 mutually cancel between 
the terms 2P2 and P 2 (the literal correspondence between 
our probabilities and those of Beenakker and Schomerus 
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is obtained by the substitution 70 — > 4a /ft). Substitut- 
ing our probabilities for this regime, we find that, in con- 
trast to the result of the Ref. [D, the leading contribution 
to 2P2 — P 2 is proportional to ry~ 2 , 

m-Pl~^T]\RT)\^) 2 . (36) 

Technically this difference arises from the fact, that 
Pi cannot compensate for the contribution from the re- 
ducible current correlator in 2P 2 , see Eq. (|2"B1 . The phys- 
ical reason for the observed distinction between our re- 
sults and the results of Ref . [l] lies in the different measure- 
ment procedure. The Glauber photodetection procedure 
implies a real counting of the photons, i.e., in addition 
to the number of photons one gains extra information 
on the time of the detection. In contrast, in our detec- 
tion scheme we do not count the plasmons — our proba- 
bilities Pi (CI) and P 2 (Cii,Ci 2 ) contain only information 
about the total number of the plasmons at the end of 
the measurement. Thus, we do not know which plasmon 
with frequency f2i or Ct 2 was detected first. As a result, 
our probability to observe two plasmons P 2 (fi) is two 
times larger than that obtained via the Glauber detection 
scheme, and hence no compensation occurs in our case. 
Instead, it is the quantity P2 = P2—P1 which exhibits the 
proper cancellation and provides a measure for the non- 
trivial correlations in the plasmon statistics in our analy- 
sis. The same situations occurs in the regime where plas- 
mons are detected in the frequency band [luq — ACl,u)o], 
ACt < o->o/2; our probability P 2 (Cl) is again twice larger 
than the result in Ref. [2} 

IV. MEASUREMENT OF THE FOURTH 
ORDER CURRENT CORRELATORS 

As mentioned in the introduction, the measurement 
of the plasmon statistics also reveals information on the 
current fluctuations in the QPC. In particular, the proba- 
bilities Pi (ft) and P 2 (Cti,Cl 2 ) provide information on the 
second- and fourth-order noise correlators at finite fre- 
quencies. Furthermore, the measurement of high-order 
current correlators, or alternatively, high-order transmit- 
ted charge cumulants, is a non trivial issue. It turns 
out, that the measurement of the irreducible correlator 
Pi provides direct access to the fourth-order charge cor- 
relator. 

The charge statistics is conveniently analyzed in a 
Gedanken experiment, where the transmitted charge is 
measured with the help of a spin- 1/2 counter^i^. Us- 
ing the expression Q t = J I(t)cLt for the transmitted 
charge, the corresponding generating function involves 
the specific time ordering 

X (X) = (21 [exp(aQ t /2)]T + [exp(iA&/2)]), (37) 

where T_ and T+ denote backward and forward time- 
ordering operators. The resulting statistics turns out to 



be binomial with the proper electron charge describing 
the transmitted carriers^. The fourth-order (zero fre- 
quency) charge cumulant is given by a weighted com- 
bination of various time orderings, 

((Qt)) = Jfif dTldT 2 dT 3 dT 4 

+4((T_(/ Tl / T2 7 T3 )J T4 ))+6«T_( )T + (I T3 J T4 ))) 

+4«j n r + (444))) + «r + (4444)))] , (38) 

and the result assumes the form 

((Qt)) = 2e 4 T(l - T)(6T 2 - 6T + 1)^11, (39) 

with the poissonian limit restored in the tunneling 
regime. 

While it is clear that this (binomial) result manifests 
itself when the charge transport is analyzed with a spin- 
1/2 counter, the question can be posed whether a more 
realistic experiment, e.g., the present plasmon detection 
experiment, can be used to measure this result. 

Indeed, the measurement of the probability correla- 
tor P 2 (Cti,Ct 2 ) tests the time-ordered fourth-order cur- 
rent correlator, reminding about the spin-1/2 detection 
scheme of Ref. [f| cf. Eqs. (J7]) and (j3"8]) . Assuming an ex- 
tended measurement where the emitted radiation is col- 
lected from a region near the QPC with a size L larger 
than the characteristic length A, the explicit calcula- 
tion of the frequency integrals in Eq. (TIT)) gives a result 
for p2(fii,fi 2 ) at low frequencies Cti,Ct 2 <C eV/H and 
|Oj — O2 j ^> 77, see Eq. (fM|). which coincides with the 
fourth-order charge correlator in Eq. (|39|) . 

P 2 (fti,ft 2 ) =a 2 T(l-T)(6T 2 -6T+l)^^i^. (40) 

hr\ \ti\ti 

We thus conclude that for an extended measure- 
ment scheme with L 3> A, the probability correla- 
tor P 2 (Cli,Cl 2 ) at low frequencies is proportional to the 
fourth-order charge cumulant ((Qt)) with t = r; _1 . 

V. CONCLUSION 

We have presented a perturbative calculation of the 
statistics of plasmon emission from electrons scattered at 
a quantum point contact. In our analysis, we determine 
the probability densities p n to find n plasmons with pre- 
scribed frequencies Cli, . . . , Cl n during a measuring time 
to ~ 1/?/; our perturbative calculation includes terms up 
to fourth order in the interaction Hamiltonian and al- 
lows us to calculate one- (pi) and two- (p 2 ) plasmon pro- 
cesses. These probability densities are related to second- 
and fourth- order current correlators and hence measur- 
ing the plasmon statistics provides also information on 
higher order current correlators. 

Our central quantity is the ireducible probability cor- 
relator P 2 (Cti,Ct 2 ) = P 2 (Oi,Q 2 ) - Pi(Oi)Pi(fi 2 ), which 
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we find to provide the most valuable information if mea- 
sured at different frequencies |Oj — £l 2 \ ^S> r\. Its sign 
provides information on the statistics of plasmon emis- 
sion, which arises from a competition between bunch- 
ing due to for multi-plasmon emission from one elec- 
tron and anti-bunching imprinted onto the plasmons 
by the regular stream of incident electrons. The char- 
acter of the emitted radiation can be tuned between 
bunched and anti-bunched by changing the voltage V 
and the transmission T across the QPC. Measuring the 
power cross-correlator in different transmission lines as 
done in recent experiment a 11 ' 12 provides experimental ac- 
cess to this quantity. Within the usual photo detection 
scheme insteadi— , the role of P% = P 2 — Py is played 
by the deviation of the variance from the Poisson value, 
Var(iV) — (N) = 2P 2 — Pf. The discrepancy in the factor 
2 in front of P2 is a consequence of the different measure- 
ment techniques, providing more detailed information in 
the photo detection scheme. This rather innocent look- 
ing difference in fact requires the definition of a different 
measurement quantity for the two cases of 'counting' (set 
of single photon measurements) and 'collecting' (single 
projection of a plasmon number state at the end); the 
former requires the analysis of Var(iV) — (iV), while the 
latter forces one to discuss P 2 . Finally, we have shown 
that the irreducible correlator ^2(^1,^2) coincides (up 
to a scale factor) with the fourth-order charge cumulant 
((Qt)) an( i nence provides practical access to this quan- 



tity. 
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APPENDIX A: FOURTH-ORDER CURRENT 
CORRELATOR 

Using the scattering matrix approach, we first calcu- 
late the irreducible fourth-order current correlator in the 
time representation, 

c(x,t) = ((i(x 1 ,t 1 )i(x 2 ,t 2 )i(x 3 ,t 3 )i(x 4 ,t 4 ))), (Ai) 

and then determine the expression (fli, fl 2 ) in 
Eq. {ISD. 

We assume a scattering process described by energy 
independent scattering amplitudes r, r, and t within a 
region near the QPC located at the origin x — 0; here, r 
and f denote reflection amplitudes for electrons coming 
from the left and right reservoirs, and t is the transmis- 
sion amplitude. Linearizing the energy-momentum dis- 
persion relation near the Fermi level, the electron current 
operator takes the form, 



I(x >0,t) = I jdede'(Ta%> + rt* a% +r*tb%, + Rbtb e ,)e z(e - e ' ){t ~^ )/ri ~ btk,e l{e - e ' ){t+ ^ )/h , (A2) 
i(x <0,t) = ~ [dede'(Ra\a e , +f*ta% + ft* b\a t , + T ^6 e ,)e j(e - e ' )(t+ ^ )/ri + at^g^-^t*-^)/^ (A3) 

I 



where a e and b e are annihilation operators for electron 
scattering states incoming from the left and right reser- 
voirs, respectively (we assume spinless electrons; T = \t\ 2 
and R = 1 — T are transmission and reflection probabil- 
ities, v F is the Fermi velocity). 

The current operator I(x, t) can be written as a sum 
I(x,t) — I~ (£ _ ) + / + (£ + ) of outgoing and incoming cur- 
rents and which depend only on the re- 
tarded variables ^ = t ± |x|/u F . Below, we concentrate 
on the current fluctuations to the right of the barrier. 
Introducing an additional index a = ± denoting the in- 
coming and outgoing currents, we rewrite the current op- 
erators in a compact form, 

= I f dedeV( £ - e ')e a A^cJ( e )Ag.c^ e '), (A4) 

ij 

where we have defined ci(e) = a e , 62(e) = 6 e , and the 



2x2 matrices A ± 

(A5) 

defined to the right of the barrier. In order to calculate 
the fourth-order current correlator C(x,t ), we have to 
average over all possible products of four current opera- 
tors, 

C S (() = ((I ai (^I a2 (^)I a3 (^ 3 )I ai (^))), (A6) 

with C(x,t) = J2s )■ Below, we use the shorthand 
Ct ' = 6 and A°" = Ai and put H = 1. Using Wick's the- 
orem and taking averages over the reservoirs, we arrive 
at the expression 
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4 i- 4 

C 5 (f) = / ]J A^e^-^e" 2 ^ 2 -^ 3 ^^^ 

^ n ' J i=i 

+ e jei( « 1 -« 4) e le2( « 2 - ?l) e le3( « 3 -« 2) e le4(?4 ^ 3) Tr{iV(e 1 )A 1 ^(e 2 )A 2 iV(e3)A3iV(e4)A4} 
_ e ^i(Ci-?4) e «2(6-«3) e «3(?3-6) e «4(C4-6) Tr {^( ei )^ l7 v( e3 )yl 3 iV(e 2 )A2iV(e4)A4} 
_ e ^i(Ci-?3) e «2(6-«i) e «3(?3-?4) e i e 4(?4-6) Tr { iV ( ei )^ 1 iV(e 2 )A 2 7V(e4)A4iV(e3)A3} 
_ e ^i(Ci-?3) e «2(6-?4) e « 3 fe-?2) e i e 4(?4-?i) Tr { iV ( ei )^ 1 iV(e 4 )A47V(e 2 )yl 2 iV(e3)A3} 
- e ^(«iH;2) e ie 2 (&-<;4) e ie3(£ 3 -^^ (A7) 

with the 2x2 matrices N(e) and N(e) defined as 

N ^=( n t ] nn(e))> m = t-N(e), (A8) 

and rii(e) and nji(e) the Fermi distribution functions of the left and right electronic reservoirs. Next, we integrate 
over the energies in Eq. (|A7[) using 



de(l-n L/R (e))e^ = . / den £/fl e)e ,gg = , — -r, A9 

' sinn|7T0(£ + zo' ) J J ' smh[7r0(£ — io")\ 

where 9 is the temperature of the fermionic reservoirs and loq — eV/ h is the voltage frequency defined by the bias 
voltage V applied to the QPC. 5' ', 5" > are regularization parameters; for an energy independent transparency, they 
are of order §' ~ hj {E c — E F ) and 5" ~ h/E F , where E c is the energy width of the conduction band and E F the Fermi 
energy. Below, we define a single regularization parameter 5 = max{(5', 5"} and assume the zero temperature limit 
9 = 0. After integration and using these simplifications we obtain the expression 



TrMSi ~ &)Aipfa - £i)A 2 p(& - &)A 3P (U - tj 3 )A 4 } + c.c. 

(2tt)4 V " (a - U - iS)(& - Ci + - 6 + - 6 + «) 

Tr{p(6 - U)Ai P (& - ~ 6)^2P(6 ~ 6)^4} + c.c. 

(6 - & - ) (6 - & + M) (6 - & + W) (& - 6 + iS) 
Tr{p(6 - - £i)A 2 pfa - &)Atp(& - £4)^3} + c.c. 

(£1 - & - «) (6 - 6 + (U - 6 + i*) & - 6 + 



(A10) 



where is the 2x2 diagonal matrix, 

p(0 = [ o "J- (AH) 

Next, we consider the irreducible contribution to the 
probability P 2 (fii,f2 2 ), see Eq. (fTTj) . This contribution 
is proportional to the frequency integral over the fourth- 
order noise correlator {uj\, ill +il 2 , w 2 ) with a specific 
kernel. Assuming a stationary situation (i.e., only rela- 
tive times are relevant) and changing from frequency to 
time variables in Eq. (TT8|) . we have to calculate the ex- 
pression 

00 

Q (4 Hn 1 ,n 2 )= J d n J d T2 c{T U n 1 +n 2 ,T 2 ) 

-00 

(• e -in 1 (r 1 +T 2 ) + e -Kn 1 r 1 +f2 2 r 2 ) +Qi ^ ( A12 j 
where we have defined 

G(n,n,7a)= [ drC(n,T,r 2 )e- mT . (A13) 



Here, C(ti,t 2 ,t 3 ) is the coordinate averaged correlator 
C{x,t) expressed in terms of the relative time variables, 



C(n,T 2 ,T 3 ) = Jd 4 xF(x) (A14) 
x C ( x ; r 1 + t 2 + r 3 , r 2 + t 3 , r 3 , ) , 



with the kernel F(x) = Ili /(^i) describing the spatial 
extension of the plasmon. 

Next, we find the non-vanishing contributions 
C s (Ti,fl > 0,t 2 ) to the correlator C(ri,fi,T 2 ) = 
^2 S (7 Q (ri, f2, r 2 ) defined in Eq. ()A13|) : these are iden- 
tified as those with a G {( ), (H ), ( h 

- -). (" - +-), ( +)> (- + +-), (+ - +-), (- + 

— +), (H h)}. The corresponding correlators can be 

written in the form, 
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(zi,z 2 )exp^ 



if2(ai|xi|+a 2 |x 2 |-a3|x 3 |-a4|z4|) ' 
2v F 



where we have introduced the new retarded variables, 

ai\xi\ - a 2 \x 2 \ 



zi(ai,a 2 ) = n 
and the functions Iq(z±, z 2 ) have the form, 



z 2 (a 3 , a 4 ) = t 2 



a3|a:3| - a4,\X4,\ 



(A15) 



(A16) 



2RT(T-R) 2 e(Lu -n)(cos[ujQ(z 1 -z 2 )/2 



sin[(cj - &)(zi+z 2 )/2] 



(z 1 + z 2 )z 1 z 2 



Z-2 



-z 2 



+2{RT) 2 Q{uj -VL/2) 



sin[(tJo - Q/2)(z! + z 2 )] 
(^i + 2: 2 )ziZ 2 



^2 



"*2 , 



( Re -i» z 2 /2 + Te i WoZ2 /2^ M^o - n)(zi + z 2 )/2] + 



(zi + z 2 )(zi-i(5)z 2 



-22 



/"++- 



7" +_+ ~ 



+2R C os[^z 2 /2] ( M(uo-n)(*i-*2)/2] 
V \Z\-Z 2 )Z\Z 2 

sin[(u/ - Jl)(zi + z 2 )/2] 



= RTQ{uj -n)e lulo{zi+Z2)/2 
= RTO{uj -n)e l " o{z2 - Zl)/2 



(z 1 + z 2 )z 1 z 2 

sin[(w - Q){z\+z 2 )/2] 
{z 1 + z 2 )z 1 z 2 

I 



-z 2 



z 2 -> -z 2 - 



-22 



-^2 



sin[(u - »)(2! + z 2 )/2] 
(z 1 + z 2 )(z 1 -iS)(z 2 -iS) 

sin[(w - Q)(zi - z 2 )/2] 
(z 1 -z 2 )(zi-i5)(z 2 -iS) 



The remaining functions 1^ can be expressed through the and t 2 in Eq. (IA12|) :. we regularize the divergent denom- 



four above, 



r-+-+ 

T+—+ 



(z 1 ,z 2 ,S) 


= rt~ 




22,- 


5), 


(zx,z 2 ,S) 


= 4- 


'(-Z 2 , 


Zl,~ 


S), 


(z 1 ,z 2 ,S) 


= 4" 


~(z 2) zi,5), 




(zi,z 2 ,5) 


= 


~{-z\, 


-z 2 , 


S), 


(zi,z 2 ,5) 






~z 2 , 


-S). 



(A17) 



inators in I^(zi,z 2 ) using the Sokhotsky formula, 



lim — !— =V-Tin$'(z), (A18) 

<5^0+ 2 ± 10 Z 



where the (5-function has to be understood as irS'(z) — 
S/(z 2 + S 2 ) with a finite width 5. The correlator 
Finally, we have to perform the time integrals over t\ Q^(f2i, fl 2 ) then takes the form 



+ Z_^( e 1 ^(l- 1 |+l^l) e -^(l a;2 l-l^l) 5 (| Xl | + | ;E2 |) + e - l ^ (|a;i|+|K4|) e-^ (|a:2| - |a:3|) lff (|x3| + |x4|) + fii «-► fi 2 ) 
i 5 (|a ;i | + |^ 2 |) 5 (|x3| + |^4|)(e-^ (|2;i| - |:E4|) e i SF (|:C2| " |x3|) +£-^(1-11+1^1)^^(1-2 1+1^1) + ^ x ^ | 

(|xi|-|x4|) e — t^-(|x 2 |-|K3|) (A19) 
I 



+ 



-2 — {RTye{uj - max{fJi,fi 2 })(w - max{Oi,n 2 }) / d 4 xF{x)e 
2-k I 



where g(x) = 1 — 0,5(x) + 0a(— x), with 0^(x) = J x dyS'(y) a Heaviside-like function with a finite width 
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X = v F S defined by the regularization parameter 5; with 
the parameters 6' ~ h/{E c — E v ) and 6" ~ h/E F) we 
obtain A ~ A F of order of the Fermi wave length. Then 
<7(|x|) = 1 for | a; | -C A and g vanishes in the opposite 
case. 

For a large collection area with A < A < L, we can 
drop all terms containing a factor g(x) in Eq. (I A19|l : the 
integration over x generates a factor ~ (v F /ft) and we 



arrive at the result Eq. (|23|) (we approximate the factors 
1/(1 + A 2 /L 2 )«1). 

Taking into account the spin 1/2 of the electron, the 
above expression has to be multiplied by a factor 2. Sim- 
ilarly, in a multi-channel situation, we have to sum the 
correlators Q$ (Q±, CI2) for all channels n with trans- 
parency T n . 
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